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We have measured time-resolved the photon storage time and the quality factor of an ultrafast photonic 
cavity using an autocorrelator. The cavity consists of a A.-thick GaAs layer sandwiched between GaAs/AlAs 
distributed Bragg reflectors and resonates at 985 nm wavelength. The inverse relative linewidth measured with 
white light reflectivity is 830, while the quality factor obtained from the time resolved measurements is 1500. 
The photon storage time in the cavity is 0.78 ps. We show that the difference between the quality factor and the 
inverse relative linewidth results from inhomogeneous broadening of the microcavity resonance due to a spatial 
gradient in the cavity layer 

PACS numbers: 



INTRODUCTION 

Since light is extremely elusive there is a great interest to 
store photons in a small volume for a certain time. Storage 
of photons in an applicable way can be achieved using solid 
state cavities. Tanabe et at. used cavities to create large pulse 
delays with small group velocities by storing light in a cav- 
ity inside a 2D photonic crystal slab fT T\. Another appli- 
cation where storage of light in a cavity plays a crucial role 
is changing the color of light as was studied by Preble et at. 
|[3l. Ultimately, with a microcavity the strong coupling regime 
of cavity quantum electro dynamics can be entered I4]|5l. In 
the strong coupling regime a cavity and a two level system 
together form a new set of states. Normal-mode splitting of 
a coupled exciton-photon mode was observed in a planar mi- 
crocavity |6|. Other interesting experiments have been per- 
formed on planar cavities, e.g Bose-Einstein condensation of 
exciton polaritons [7] and the investigation of the limitations 
of a scanning Fabry- Perot interferometer 

An important characteristic parameter of a cavity resonance 
is the storage time of light Zcav The storage time is defined 
by the response of the cavity resonance to a Dirac pulse. Ex- 
citation of the electromagnetic field in a cavity was studied in 
1^. The response to the Dirac pulse is given by an exponential 
decay of the intensity I(t) in the cavity resonance ifTOl : 

/(f) = v-'/''■''^ (1) 

with Iq the initial intensity that the pulse stores in the cavity. 
However, in more complex cavities the behavior of the cavity 



can be very different from the single exponential case ifTTl . To 
compare cavities independent of their resonance frequencies 
COo, the widely used figure of merit is the resonance quality 
factor Q, which is defined as: 

Q = T,.„vCOo. (2) 

Physically, the quality factor is proportional to the ratio be- 
tween the total energy stored and the energy lost per cycle. At 
optical frequencies a cavity with a feasible high quality factor 
of Q= 10^ is relatively slow with a response time in the or- 
der of nanoseconds. A cavity with a moderate quality factor 
Q — 1000, however, is fast with a response time of picosec- 
onds. The picosecond timescale allows ultrafast access and 
storage of light in the cavities. 

A common procedure to estimate the quality factor of a cav- 
ity is to measure a transmission or reflectivity spectrum and 
extract Q from the relative linewidth of the cavity resonance 
||3]|4l|5l[l2l[l3][l4l[l5J. For a single resonance without de- 
phasing, one can use the Wiener-Khintchine theorem, which 
relates the field autocorrelate to the intensity spectrum, to ob- 
tain 



However, if there is significant dephasing, e.g. due to inho- 
mogeneous broadening or thermal noise, Aco will in general 
be larger and Q > coo/Aco. 

From many resonating systems in condensed matter and 
solid state physics, it is known that besides homogeneous 



2 



broadening there is also the possibility of inhomogeneous 
broadening of a resonance llT6l|T7| . In the case of an ensem- 
ble of resonators inhomogeneous broadening of a resonance 
results from inhomogeneities in the resonance frequency. If 
the resonance frequency is different for each resonator the 
linewidth of the ensemble is broader than the linewidth of a 
single resonator and the ensemble linewidth is typically deter- 
mined by the distribution of resonance frequencies. In the case 
of inhomogeneous broadening the linewidth will only give a 
lower boundary for the range of possible Q values. The true 
quality factor must in this case be determined from dynamic 
measurements. 

A dynamic measurement to determine the quality factor is 
a cavity ring down experiment as was treated in 1 18 1. In this 
case a cavity is excited by a pulse and the intensity emitted 
from the cavity is measured as a function of time. In the case 
of storage times in the order of nanoseconds and very high 
quality factors {Q = 10*') time correlated single photon count- 
ing can be used to determine the storage time 1 1 1. In our case 
of ultrafast cavities that decay on a ps timescale with mod- 
erate quality factor {Q = 1000), an intensity autocorrelation 
function is the method of choice for determining the quality 
factor. 

The normalized correlation function that is measured is an 
intensity autocorrelation function G^, which is given by ||T9l 

G\x) = ^ ^ , '—, (4) 

^0 

where T is the delay time between the pulses from each of 
the interferometer branches, /q is equal to maximum value of 
the unnormalized autocorrelation value, and /(f) is the time 
dependent intensity. There is no phase in equation |4j which 
means that this is the proper autocorrelation function, also 
in case of dephasing. The autocorrelate has its maximum 
at delay T = 0, when the pulses in the two branches of the 
Michelson interferometer overlap. For example the autocor- 
relate of a Gaussian pulse is given by a Gaussian shape, where 
the width of the input pulse Zip and the autocorrelate are re- 
lated as lac — \/i^)'^ip- From the autocorrelate of a pulse 
stored in the cavity resonance, the storage time can be found 
from the full width at half maximum ZfWHM of G^, with 
"^cav ~ 0.63zfwhm- 

EXPERIMENTAL 

Our structure is a planar cavity that consists of a GaAs A,- 
thick layer (277 nm thick), sandwiched between two Bragg 
stacks. One Bragg stacks consists of 12 and the other Bragg 
stack consists of 16 pairs of A,/4-thick layers of nominally 
pure GaAs or AlAs. The same structure was studied in Ref. 
[151. The sample was grown at CEA in Grenoble by means of 
molecular beam epitaxy at 550"C 1201 . For experiments out- 
side the present scope the sample was doped with lO'^cm^^ 
InGaAs/GaAs quantum dots, which hardly influence our ex- 



periment ll28l . There is a spatial gradient in the cavity thick- 
ness of ^ = 5.64 nm/mm 121]. The spatial gradient results 
in a position dependent resonance frequency. In our measure- 
ments we average the transmitted intensity over the area of 
the focal spot. The different resonance frequencies cause the 
resonance to broaden inhomogeneously. 

White-light reflectivity and transmission were measured 
with a broadband white-light spectrometer setup with a spec- 
tral resolution of about AA- = 0.2 nm ll22l . The transmission 
spectrum was measured with a collimated beam with a diam- 
eter of 2 mm. The reflectivity spectrum was measured with 
a glass objective with a numerical aperture NA = 0.05 and a 
focus diameter of 100 fjin. The reflectance spectrum of a gold 
mirror was used as a reference. 

For pulse transmission and the intensity autocorrelate, we 
used a Titanium Sapphire laser that emits Zip = 0.115 ps 
pulses at A, = 800 nm at a repetition rate of 1 kHz (Hurri- 
cane, Spectra Physics). The laser drives an optical parametric 
amplifier (OPA, Topas 800-fs, Light Conversion), which gen- 
erates the pulses used to probe the photonic cavity. The center 
wavelength of the OPA pulses can be tuned between 450 nm 
and 2400 nm. We used a fiber optic spectrometer (USB2000, 
Ocean Optics) to measure transmission spectra of the fem- 
tosecond pulses. We measured with an unfocused collimated 
beam with a spot diameter of 2 mm, and a numerical aperture 
NA = 10^"*. The intensity autocorrelation function was mea- 
sured using a Pulse Check autocorrelator (APE GmbH). The 
autocorrelator consists of a Michelson interferometer with a 
scanned delay path and a nonlinear crystal that generates sec- 
ond harmonic light. The autocorrelator has a maximum range 
of 15 ps with a resolution of 1 fs. We used the same beam 
parameters as in transmission. The intensity on the sample is 
100 kWcm^^, sufficiently low to avoid non-linear effects. 

Simulations were performed with the finite-difference time- 
domain (FDTD) method using a freely available software 
package with subpixel smoothing for increased accuracy [231 . 

EXPERIMENTAL RESULTS 

In figure[T]A we show the reflection and transmission spec- 
tra of the planar cavity. A prominent stopband with a reflec- 
tion of 100 % and a transmission of % is visible. Outside the 
stopband a Fabry-Perot fringe pattern is visible, while inside 
the stopband a narrow trough in reflection and a narrow peak 
in transmission mark the position of the cavity resonance. An 
effect of the spatial gradient in the cavity thickness is visible 
in the spectra in figure [T]B: The frequencies of the peak and 
trough, which are measured at different sample position, dif- 
fer slightly. Reflectivity and transmission measurements on 
the same spot are shown as an inset in figure[T]B. The trough 
and the peak are clearly at the same wavelength as expected. 

The solid line in figures [T] A and[T]B represents a transfer 
matrix (TM) calculation, with fixed complex input parameters 
nooAs Ii24i and n^iAs |25|. The thickness of the A/4 layers 
(dcaAs = 70.2 nm and cIaias — 83.2 nm) and the thickness of 
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FIG. 1: (A) Linear reflectivity and transmission spectrum of the 
GaAs/AlAs microcavity. The solid line represents the fit with a Trans- 
fer Matrix (TM) model. A stop band is apparent in both reflection and 
transmission; the trough in the reflectivity spectrum and the peak 
in the transmission spectrum reveal the presence of the cavity. (B) 
Zoom-in of (A). From the linewidth of the trough and peak, an in- 
verse relative linewidth of 830 was found in reflection and transmis- 
sion. The resonance is slightly shifted between the transmission and 
reflection measurement due to realignment of the sample between the 
measurements. The inset in B shows a reflection and a transmission 
spectrum, measured at the same position. Here the trough and peak 
are clearly at the same wavelength. 



the cavity (dcm = 277 nm) were obtained by fitting the resuhs 
of the calculations to the measured spectrum. These values 
are in agreement with expected values from the fabrication 
process. The calculation fits well with respect to frequency 
and amplitude. The reflectivity of the measured stopband is 
higher than the calculated value of 100 % because of a small 
systematic error in the gold reference spectrum. 

It is apparent from figure [T]B that the calculated linewidth 
of the cavity resonance is narrower than the measured 
linewidth. We attribute this discrepancy to inhomogeneous 
broadening of the measured linewidth, due to the spatial gra- 



dient in the cavity layer thickness. With the 100 ijm diameter 
spot we average over different positions and therefore over 
different resonance frequencies. Broadening due to a spread 
in wavevectors can be neglected since the numerical aperture 
of the impinging beam was made very small {NA < 0.05), as 
opposed to 1 15 1, where a high NA was used. We find that the 
relative linewidth in both reflection and transmission equals 
^ = 830, with AX the full width at half maximum (FWHM) 
and A-o the resonance wavelength. The transfer matrix cal- 
culation yields an inverse relative linewidth of 1640 ± 100, 
about double the value of the inverse linewidth measured with 
white-light spectroscopy. 

We measured the intensity autocorrelation traces to deter- 
mine the true storage time and Q of the cavity resonance with 
a time-resolved measurement. Figure [2] shows the autocorre- 
lation traces at values of the center wavelength of the OPA, 
XopA = 930 nm (A), Xqpa ^ 985 nm (B) and Xqpa = 1070 
nm (C). All figures show that the pulses that are transmitted 
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FIG. 2: Normalized intensity autocorrelation traces of pulses trans- 
mitted through a planar cavity at different OPA wavelength settings: 
930 nm (A), 985 nm (B) and 1070 nm (C). The autocorrelation traces 
of the input pulses are given by the circles, while the autocorrelation 
traces of pulses transmitted through the cavity are offset by 0. 9 and 
given by squares. The dashed and solid lines areflts to the autocorre- 
lation traces, without and with sample respectively. The shape of the 
autocorrelation trace is Gaussian for the pulses from the OPA. The 
pulses that are on resonance with the cavity show an autocorrelate 
that agrees very well with the autocorrelation trace from the damped 
oscillator model (B). The shape of the pulses transmitted through a 
non-photonic range of the sample remains Gaussian. 

through the sample are broader than the input pulses. The 
width of the input pulses is T,p = 0. 1 15 ps and the shape Gaus- 
sian, which we expect from the specifications of our laser sys- 
tem. The transmitted pulses are broadened by dispersion in 
the off resonance cases (A) and (C). In the case of figure|2]B 
the broadening is the result of the storage of the photons in the 
cavity. 

The shape of the autocorrelation trace of the transmitted 
pulses is Gaussian for pulses transmitted outside the stop 
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band, as expected. The autocorrelation traces measured on 
resonance with the cavity (B) are non-Gaussian. This is typ- 
ical for autocorrelation traces near the cavity resonance, be- 
cause of the exponential decay of the energy stored in the 
cavity. The autocorrelation traces calculated with a damped 
oscillator model is shown in figure (B) and fits the experimen- 
tal data very well. From the width of the autocorrelation trace 
on resonance {Ifwhm = 1.1 ps), we conclude that the true 
storage time of our cavity is icav = 0.78 ± 0.05 ps. 

To farther analyze the autocorrelation traces we plot the 
full width at half maximum of the measured autocorrelation 
traces. The results are shown in figure |3] A as a function of 
center wavelength of the laser. The width of the autocorrelate 
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FIG. 3: (A) Pulse width as a function of wavelength setting. The 
pulses that passed through the sample (squares) are broadened with 
respect to the input pulses (circles). The solid/dashed line represents 
the FWHM of a 15 nm spectrally wide input pulse that is transmitted 
through a 350 pm GaAs wafer that is used as a substrate. (B) Nor- 
malized spectrum transmitted through the microcavity, and reference 
spectra of pulses directly from the OPA. For 935 nm wavelength set- 
ting (green squares) the transmitted spectrum consists of the tail of 
the input spectrum that continues in the region on the blue side of the 
stop band (blue circles). 



of pulses without the sample is Xac = 0.115 x ^/2 ps and es- 
sentially independent of laser wavelength, as expected from 



the OPA specifications. In the presence of the sample, we ob- 
serve a more complex dependency on the wavelength, with 
three regimes: Transmission on resonance, transmission out- 
side the stopband, transmission inside the stopband. Near the 
cavity resonance the width of the autocorrelate increases dras- 
tically to Ifwhm — 1-1 ps. The width of the autocorrelate at 
the cavity resonance is attributed to the storage of light in the 
cavity: The storage time of the cavity icav — 0.78 ± 0.05 ps 
and the quality factor is equal to 1500 ± 100. 

Outside the stopband the pulses are broadened. The width 
of the autocorrelate outside the stop band is about 0.2 x \/2 ps. 
We attribute the broadening outside the stopband region to dis- 
persion in the GaAs substrate. From figure [3}\ it can be seen 
that the width of the autocorrelation traces matches well the 
expected width for a pulse transmitted through a GaAs wafer 
||26| . The expected width is calculated for a GaAs wafer with 
a thickness of 350 /jm, from the dispersion given by Blake- 
more Il24l . 

In figure |3] A, we observe datapoints inside the stopband, 
where a transmission of % is expected. We measure values 
for the width that are close to the values outside the stopband. 
The situation in this case is sketched in figure [3] B where we 
see the transmitted spectrum with and without sample. We 
observe that the blue part of the spectrum is transmitted, which 
means that the measured with of the intensity autocorrelate is 
the value for the blue side of the stopband. 



MODELING 

To obtain a physical picture of the decay mechanism inside 
the cavity and to verify what the true Q is, we model the be- 
havior with a damped harmonic oscillator We furthermore 
performed FDTD calculations to calculate the Q of the cav- 
ity in the ideal case and to check the validity of the harmonic 
oscillator model 129|. 

The response of a damped harmonic oscillator with Q = 
1450 to a Gaussian input pulse with a width of T,p = 0.2 ps 
is shown in figure |4] together with the Gaussian input pulse 
and the cavity response as calculated with FDTD. No disper- 
sion and no absorption was taken into account for the FDTD 
calculations. In the harmonic oscillator case and in the FDTD 
case the intensity decays exponentially and with the same rate. 
Therefore, we conclude that the harmonic oscillator is a suit- 
able model to describe in a simple way the decay of the mi- 
crocavity. Furthermore the quality factor of the cavity without 
absorption and dispersion is equal to Q= 1450 ± 100. 

With the damped harmonic oscillator model we have cal- 
culated the autocorrelates that are shown together with the 
measured data in figure |5] for a quality factor of the damped 
harmonic oscillator Q = 1500 and Gaussian input pulse with 
width 0.12 ps. Figure |5] shows a very good agreement be- 
tween the measured autocorrelation trace and the calculated 
autocorrelation trace for Q = 1500. The autocorrelate obtained 
from this model has a FWHM of 1 .2 ps, which is in very good 
agreement with the measured value of 1.1 ps. 
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FIG. 4: Response of a damped harmonic oscillator with a quality 
factor Q = 1450 (dashed) to a Gaussian input pulse (solid). The 
output pulse is the input pulse convoluted with the impulse response 
of the oscillator. The symbols represent the resulting decay of the 
intensity in the cavity, as obtained from the FDTD calculation. 
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FIG. 5: Autocorrelate of pulse transmitted at the cavity resonance 
(circles). The solid line represents the result from a simple oscillator 
model with Q = 1500. The input pulse duration is 0.12 ps. 



TABLE I: Overview of inverse linewidtlrs and cavity lifetimes mea- 
sured and calculated with presented methods. 



Method 



AX 



Measurement 



Autocorrelator 
Refl./Trans. 



1500 ±100 
830 ±50 



0.78±0.05 
0.43 ±0.02 



Theory 



Transfer Matrix 
FDTD 



1640 ±100 
1450 ±100 



0.93±0.05 
0.75 ±0.05 



measurements. The difference results from inhomogeneous 
broadening due to the spatial gradient in the thickness of the 
cavity layer. For a quality factor of 1500 we expect a width 
of 0.64 nm, while we measure a width of 1.2 ±0.1 nm. Be- 
cause of the focus diameter of 100 /jm and the spatial gradient 
of 5.64 nm/mm, we expect a broadening of 0.56 nm. We find 
a total width of 1 .2 ± 0. 1 nm if we add the broadening to the 
unbroadened width. The total width of 1 .2 ± 0. 1 nm is in per- 
fect agreement with the measured value of 1 .2 ± 0. 1 nm. The 
FDTD calculation agrees very well with the measured value 
and the value obtained from the transfer matrix calculation. 
No absorption is taken into account in the FDTD calculation, 
which is the case for the TM model. 

The inverse relative linewidth of 830 that we find from the 
transmission and reflectivity measurements is the result of in- 
homogeneous broadening of the resonance. The origin of this 
inhomogeneous broadening is most likely the spatial gradi- 
ent in the A,-thick GaAs layer. Because there is inhomoge- 
neous broadening, the planar microcavity should be viewed 
as a static ensemble of microcavities of which each resonance 
frequency is slightly shifted L16J . This is in agreement with 
the results presented in lIZTl . where the spatial extent of modes 
was investigated. The effective radius of the mode is given by 

'"e// ~ (G/^'^)!^/")^- °™ ''^^^ f'li'l ^ spatial extent in 
the order of 5 /jot, which is much smaller than the diameter of 
the probe beam. In general our results show that the true qual- 
ity factor of a planar microcavity indeed can only be obtained 
from a time-resolved measurement. 



CONCLUSION 



Table [ijpresents the results of the various methods to deter- 
mine the quality factor presented in this work. We see that 
the value for the inverse linewidth found from the transfer 
matrix model agrees well with the value found from the au- 
tocorrelator measurement. The slight discrepancy might be 
attributed to minor irregularities in the structure. We con- 
clude that the quality factor of the cavity is 1500 ± 100 and 
the storage time of light 0.78 ± 0.05 ps. The inverse linewidth 
measured in reflection and transmission with white light spec- 
troscopy is much smaller than the value from autocorrelator 



In the case of an inhomogeneously broadened resonance 
we have shown that the intensity autocorrelate can be used 
to determine the storage time of a cavity resonance. For an 
inhomogeneously broadened microcavity resonance we have 
measured both the spectral width and the intensity autocorre- 
lation trace. The intensity autocorrelation trace yields a value 
of the quality factor that agrees well with the values found 
from transfer matrix and FDTD calculations. The spectral 
width is affected by inhomogeneous broadening and leads to 
the wrong value for the quality factor. 
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